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A recent mega-analysis combining genome-wide association study data revealed that a variant of microRNA 1 37 (MIR 1 37) exhibits the 
most significant association with schizophrenia. Other biological evidence also consistently suggests that MIRI 37 may have a pivotal role 
in the pathogenesis of schizophrenia. However, the underlying neural mechanism remains unclear. As the disrupted dorsolateral 
prefrontal cortex (DLPFC) coupling with the hippocampal formation (HF) has been widely observed in schizophrenia patients, DLPFC- 
HF dysconnectivity can therefore be thought of as a pivotal intermediate phenotype that links genetic variants of psychiatric risk genes to 
schizophrenia. This study used resting-state functional magnetic resonance imaging to test whether the MIRI 37 variant (rs 1 625579) 
impacts DLPFC-HF functional connectivity and cognitive performance in 290 young, healthy Han Chinese individuals. To identify 
functional connectivity between DLPFC and HF, a seed-based functional connectivity analysis was used. The association between 
DLPFC-HF connectivity and working memory performance was further examined in individuals with different MIRI 37 genotypes. The 
individuals who are homozygous for the MIRI 37 risk allele (TT), which confers a high risk for schizophrenia, exhibited significantly 
different DLPFC-HF functional connectivity compared with TG individuals. Moreover, the DLPFC-HF connectivity could predict the 
working memory performance in MIRI 37 TG individuals, but not in TT individuals. The current findings obtained in a large sample of 
healthy participants identified potential neural mechanisms linking MIRI 37 with the risk of developing schizophrenia via the intermediate 
phenotype of DLPFC-HF connectivity. 
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INTRODUCTION 

A recent mega-analysis combining genome-wide association 
study data revealed that microRNA 137 (MIRI 37) rsl625579 
exhibits the most significant association with schizophrenia, 
and that four predicted targets of MIRI 37 also display 
genome-wide significance (Ripke et al, 2011). Several 
schizophrenia-associated genes (CSMD, C10orf26, CAC- 
NA1C, TCF4, and ZNF804A) were further experimentally 
validated as miR-137 targets (Kim et al, 2012; Kwon et al, 
2013). Furthermore, MIR137 rsl625579 were recently 
reported to be associated with impaired cognitive perfor- 
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mance of schizophrenia patients (Cummings et al, 2013; 
Green et al, 2013). Moreover, miR-137 has a role in 
regulating adult neurogenesis (Szulwach et al, 2010) and 
neuronal maturation (Smrt et al, 2010), which is consistent 
with the neurodevelopmental hypothesis for the etiology of 
schizophrenia. All of this evidence suggests that MIRI 37 
may have a pivotal role in the pathogenesis of schizo- 
phrenia. However, the neural mechanisms underlying this 
association are largely unknown. 

Convergent evidence has indicated that the dorsolateral 
prefrontal cortex (DLPFC) is one of the most important 
regions involved in the pathogenesis of schizophrenia. 
Functional magnetic resonance imaging (fMRI) studies have 
demonstrated abnormal DLPFC activation in schizophrenia 
patients and their unaffected relatives during cognitive tasks 
(Callicott et al, 2003; Glahn et al, 2005; Seidman et al, 2006; 
Karlsgodt et al, 2007; Becker et al, 2008). More intriguingly, 
a disrupted DLPFC coupling with the hippocampal forma- 
tion (HF) has been observed in schizophrenia patients 
(Meyer-Lindenberg et al, 2005), and similar findings have 
been reported in individuals at a high genetic risk for 
developing schizophrenia (Esslinger et al, 2009; Rasetti 
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et al, 2011; Paulus et al, 2013, 2014). Previous studies have 
also shown that the interaction between the hippocampus 
and the prefrontal cortex supports working memory (WM), 
long-term memory, and other cognitive abilities that are 
abnormal in schizophrenia patients (Axmacher et al, 2008; 
Cohen, 2011; Hanlon et al, 2012; Supekar et al, 2013). 
Resting- state functional connectivity, which is measured 
by resting-sate fMRI, examines temporal correlations of 
intrinsic low- frequency fluctuations in the blood oxygena- 
tion level-dependent (BOLD) signals across brain regions, 
and has been proven useful for mapping brain networks 
(Biswal et al, 2010; van den Heuvel and Hulshoff Pol, 2010; 
Rosazza and Minati, 2011) and may predict individual 
cognitive performance (Kelly et al, 2008; Hampson et al, 
2010; Supekar et al, 2013; Unschuld et al, 2013). Resting- 
state DLPFC-HF connectivity is therefore thought of as a 
pivotal intermediate phenotype linking genetic variants of 
psychiatric risk genes to schizophrenia and related cogni- 
tive performance. 

Previous experimental evidence demonstrated that miR- 
137, as a brain-enriched microRNA, has important roles 
during both the early and late stages of adult neurogenesis, 
including the proliferation and differentiation of adult 
neural stem cells, the regulation of dendritic morphogen- 
esis, and synaptic plasticity (Silber et al, 2008; Siegel et al, 
2009; Smrt et al, 2010; Szulwach et al, 2010). Based on such 
biological role of miR-137, we may speculate that miR-137 
influences the development and plasticity of hippocampal- 
cortical circuitry. Meanwhile, the dys connectivity hypoth- 
esis of schizophrenia has been widely proven in various 
neuroimaging studies (Pettersson-Yeo et al, 2011). The 
connectivity-based neuroimaging intermediate phenotype, 
acting as a bridge that well links the biological function of 
MIR 137 to the neural mechanisms of schizophrenia, may 
help reveal specific mechanisms that link MIR 137 to the risk 
of schizophrenia. Therefore, by introducing the intermedi- 
ate phenotype of DLPFC-HF connectivity (as measured by 
resting-state fMRI), we aimed to provide a novel way to 
understand the mechanisms of the association between 
MIR 137 and schizophrenia. 

Here, we hypothesized that individuals who are homo- 
zygous for the MIR137 risk allele (TT), which confers high 
risk for schizophrenia (Ripke et al, 2011), would display 
disrupted DLPFC-HF functional connectivity and different 
effects on WM performance compared with other indivi- 
duals. The present study applied an imaging genetic strategy 
in 290 genotyped healthy Han Chinese participants to 
investigate the hypothesis. 

MATERIALS and METHODS 
Participants 

We totally recruited 323 young healthy Han Chinese 
participants (157 males and 166 females, mean age = 22.7 
years, range =18-31 years, all right-handed). After they 
were provided with a complete description of the study, all 
participants provided written informed consent. The pro- 
tocol was approved by the Ethics Committee of Tianjin 
Medical University. All of the participants were carefully 
screened to ensure that they had no history of neurological 
or psychological diseases in the subjects or their first-degree 



relatives, psychiatric treatment, or drug or alcohol abuse, 
traumatic brain injury, or visible brain lesions on conven- 
tional MRIs. All subjects were examined using the Chinese 
Revised Wechsler Adult Intelligence Scale (WAIS-RC) 
(Gong, 1982). After excluding three subjects with missing 
MIR137 genotype data and 30 subjects who lacked sufficient 
fMRI data, a total of 290 subjects were included in the data 
analysis of the current study. 

MIR137 Genotyping 

We extracted genomic DNA from whole blood using the 
EZgene Blood gDNA Miniprep Kit (Biomiga, San Diego, 
CA). MIR137 rsl625579 and ZNF804A rsl344706 were then 
genotyped using the iPLEX Gold Assay (Sequenom, San 
Diego, CA), following the manufacturer's instructions and 
using primers (Supplementary Table 1), and designed using 
Sequenom Assay Design 3.1 software. The accuracy of geno- 
typing was assessed by analyzing every sample in duplicate. 
Alleles were automatically called with Sequenom's MassAR- 
RAY Typer 4.0 software and verified by two independent 
reviewers. Three subjects with missing genotype data were 
excluded from further analysis. MIR137 GG genotypes were 
not identified in our current population. The reference 
frequency of rs 1625579 in dbSNP database (http://www.ncbi. 
nlm.nih.gov/projects/SNP/snp_ref.cgi?rs= 1625579) supports 
this phenomenon. 

Image Data Acquisition 

All of the subjects were scanned in the same Signa HDx 3.0 T 
magnetic resonance scanner (GE Healthcare, Milwaukee, WI) 
and were subjected to a high-resolution 3D Tl -weighted 
brain volume (BRAVO) MRI sequence with the following 
parameters: repetition time (TR) = 8.1ms, echo time 
(TE) = 3.1 ms, flip angle = 13°, voxel size = 1x1x1 mm 3 , 
and 176 sagittal slices. Resting-state functional imaging was 
subsequently performed using a single-shot-gradientecho- 
echo planar imaging (SS-GRE-EPI) sequence that is 
sensitive to BOLD contrast with the following parameters: 
TR = 2000 ms, TE = 30 ms, field of view (FOV) = 240 x 
240 mm 2 , matrix = 64 x 64, flip angle = 90°, voxel size = 
3.75 x 3.75 x 4.0 mm 3 , 40 slices, and 180 volumes. Before 
the scanning, all subjects were instructed to move as little as 
possible, keep their eyes closed, think of nothing in 
particular, and avoid falling asleep. Then, we asked subjects 
whether they had fallen asleep during and after the scanning 
to confirm that none of them had done so. 

fMRI Preprocessing 

We first excluded 14 subjects owing to apparent interslice 
motion artifacts in their raw fMRI data, which was evaluated 
via visual inspection by two radiologists who were blinded 
to genotype (C Yu and W Qin). Further data preprocessing 
was completed in the remaining subjects using SPM8 
(http://www.fil.ion.ucl.ac.uk/spm/software/spm8/) and in- 
house software. The following preprocessing steps were 
performed: (1) discarding of the first 10 volumes; (2) slice 
timing correction; (3) realigning the volumes to the first 
volume to correct for interscan movements; (4) spatially 
normalizing to a standard EPI template; (5) spatially 
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smoothing using a 6 mm Gaussian kernel; (6) linear 
regression to remove the influence of head motion, whole 
brain signals and linear trends; (7) temporal band-pass 
filtration (0.01-0.08 Hz). Specifically, sixteen subjects who 
exhibited a maximum displacement in any of the cardinal 
directions (x, y, z) of > 2 mm, or a maximum spin (x, y, z) 
of >2° were excluded from subsequent analyses. 

Functional Connectivity Analyses 

For each individual, we calculated voxel-wise DLPFC func- 
tional connectivity maps throughout the whole brain by 
taking the right DLPFC as the region-of-interest (ROI). The 
mask of right DLPFC was defined as a spherical region with a 
radius of 6 mm at the center of the MNI coordinate (42, 33, 
33), according to a recent meta-analysis study (Rottschy et al, 
2012). Each functional connectivity map for each individual 
was computed by averaging the BOLD time series in the ROI 
and then computing the Pearson's correlation coefficient 
between the average time series and those from each voxel 
in the brain. The resulting correlations were transformed 
to approximate a Gaussian distribution using Fisher's z 
transformation. Thus, the whole-brain functional connec- 
tivity map of the right DLPFC was created for each subject. 

WM Performance 

The WM task was performed on a computer in a quiet room 
outside the MRI scanner before performing magnetic 
resonance scanning for each subject. WM was assessed 
with a letter 2-back task as described previously (Ding et al, 
2012). Briefly, each participant viewed a series of letters that 
were presented sequentially and was asked to judge whether 
the letter on the screen was identical to the one presented 
two trials back. The task consisted of three blocks (30 trials 
each). Before the experiment, participants were verbally 
instructed and had a practice block, to make sure they have 
understood the task before taking the formal tests. The 
accuracy was used as the index of WM performance in the 
current study. 

Statistical Analysis 

We used a Pearson's j 1 test to check for gender differences 
and two-sample f-tests to check for differences in age, 
education levels, and intelligence quotient (IQ) scores 
between the two genotype groups. Voxel-wise one-sample 
f-test for the DLPFC functional connectivity map were 
performed in the whole group, the TT group and TG group 
(family-wise error (FWE) corrected p<0.05). Voxel-wise 
two-sample f-tests with two covariates (age and gender) 
were implemented in SPM8 to map group differences in the 
functional connectivity of the right DLPFC between MIR137 
TT homozygotes and TG individuals. The FWE correction 
was performed using the small volume correction con- 
strained to the predefined HF ROI (corrected P<0.05), 
which was created using the anatomical labels provided 
by the Wake Forest University PickAtlas (http://www.fmri. 
wfubmc.edu/downloads). The hippocampal volume was 
determined using Freesurfer software (Athinoula A Marti- 
nos Center for Biomedical Imaging, Charlestown, MA) 
(Fischl et al, 2002) and manually edited according to the 
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Tl -weighted structural MRI data. The hippocampal volume 
was then added as a new covariate for another round of 
two-sample f-tests. Fisher's z-transformed correlation 
coefficients that were significantly different between two 
genotypes in the left HF (right DLPFC functional con- 
nectivity analyses) was used as a measure of connectivity. 
For each genotype group, Pearson's correlation analysis was 
further performed to test the association between DLPFC- 
HF functional connectivity and WM performance. 

Finally, to eliminate the possibility of false positives due 
to significantly different sample sizes in the two genotype 
groups, we also performed 10 000 iterations of simulation- 
based comparisons by randomly resampling 41 of the 249 
TT subjects and comparing the random subsample with the 
same number of TG subjects. 

RESULTS 

After excluding three subjects with missing MIR137 
genotype data and the 30 subjects without sufficient fMRI 
data from the group of 323 participants, a total of 290 
subjects were finally included in the current study. The 
distribution of MIR137 genotypes (TT: N = 249; TG: N = 41; 
GG: N = 0) did not significantly deviate from Hardy- 
Weinberg equilibrium (# 2 = 1.68, P = 0.20). No significant 
differences in age, gender, years of education, IQ scores, and 
WM performance were observed between the TT and TG 
individuals (Table 1). 

In the 290 healthy participants included in the study, we 
first investigated the influence of MIR 137 polymorphism on 
the connectivity of the DLPFC to the whole brain. When the 
DLPFC functional connectivity map of each individual was 
obtained using the right DLPFC as the seed region, we 
performed one-sample £-test for the whole group, TT group, 
and TG group, and found that the three groups exhibited 
similar DLPFC functional connectivity patterns (Figure 1). 
However, individuals homozygous for the MIR 137 risk 
allele (TT) exhibited significantly different functional con- 
nectivity between the DLPFC and the left HF (peak voxel 
MNI coordinate: x=-27, y=-12, z= - 18; T = 3.43; 
cluster size = 51) compared with heterozygotes (TG) 
(Figures 2a and b). This difference survived the FWE 
multiple comparisons correction when the voxel-wise com- 
parisons were corrected using the small volume correction 
within the predefined HF region, as applied in previous 
studies (Esslinger et al, 2009; Rasetti et al, 2011; Paulus et al, 
2014), and remained significant after controlling for indivi- 
dual hippocampal volume. Interestingly, the DLPFC-HF 



Table I Participant Demographic Data by MIR 1 37 Genotype 





TT 


TG 


P- value 


N 


249 


41 




Male/female 


1 15/134 


19/22 


0.98 


Age (years) 


22.74 (2.40) 


22.85 (2.42) 


0.78 


Years of education 


15.50 (2.63) 


1 5.70 (2.22) 


0.65 


Full scale intelligence quotient 


1 17.29 (9.25) 


1 15.71 (7.15) 


0.30 


Working memory performance 


80. 1 1 (4.85) 


78.88 (5.29) 


0.14 



Values are shown as mean (SD). 

Major and risk allele = T; minor allele = G. 
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Figure I The dorsolateral prefrontal cortex (DLPFC) functional connectivity pattern obtained by one-sample t-tests for the whole group (a), microRNA 
137 (MIR 1 37) TT homozygotes (b) and TG heterozygotes (c). Threshold was set as whole brain family-wise error (FWE) corrected P<0.05. 




DLPFC-MF Functional Connectivity 

Figure 2 Effect of microRNA 137 (MIR 1 37) rs I 625579 genotype on dorsolateral prefrontal cortex-hippocampal formation (DLPFC-HF) functional 
connectivity and working memory performance, (a) DLPFC functional connectivity analyses using the right DLPFC (seed voxel MNI coordinate: x = 42, 
y = 33, z = 33) as the seed region revealed a increased DLPFC connectivity with the left HF (peak voxel MNI coordinate: x = — TJ,y= — 1 2, z = — I 8) in 
the homozygotes for the risk allele (TT) compared with TG individuals (family-wise error (FWE) corrected p<0.05). (b) Independent region-of-interest 
(ROI) analyses by taking the structural HF as ROI further supported that DLPFC-HF functional connectivity was different between the two different 
genotype groups (mean ± SD). (c) Individual DLPFC-HF connectivity negatively correlated with individual working memory performance in TG individuals 
(right red scatter diagram), whereas the correlation did not exist in TT individuals (left blue scatter diagram). 
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connectivity was inversely correlated with WM performance 
in TG individuals, whereas such association disappeared in 
individuals homozygous for the MIR137 risk allele (TT) 
(Figure 2c). This finding demonstrated that depending on 
the MIR137 genotype, DLPFC-HF connectivity may affect 
WM performance in a totally different way. 

To investigate the influence of different sample size 
between groups on the results, we performed simulation- 
based comparisons by randomly resampling 41 of the 249 
TT subjects and comparing the DLPFC-HF connectivity 
of the random subsample with the same number of TG 
subjects. A total of 8879 out of the 10 000 comparisons made 
between groups of equal sample sizes indicated that the 
DLPFC-HF connectivity in the TT subsample was signifi- 
cantly higher than in the TG group, whereas none of the 
comparisons showed the inverse result (two-sample f-tests, 
P<0.05). 



DISCUSSION 

In our imaging genetic analysis, based on resting-state fMRI 
data obtained in a large sample of healthy Han Chinese 
subjects, we found that individuals homozygous for the 
MIR 137 risk allele (TT) display significantly different 
DLPFC-HF functional connectivity compared with other 
individuals. Moreover, higher DLPFC-HF connectivity may 
predict worse WM performance in TG individuals; however, 
the prediction capability did not exist in TT individuals. 
These findings confirmed our hypothesis that MIR137 
impacts DLPFC-HF coupling and its association with WM 
capacity. The current findings provide important implica- 
tions regarding individual differences in the genetic 
contribution of MIR 137 to brain function and the risk of 
developing schizophrenia. 

One of the major findings of the present study was the 
influence of MIR137 on DLPFC-HF connectivity. As MIR137 
rs 1625579 was first reported to exhibit a strong association 
with schizophrenia (Ripke et al, 2011), MIR137 has been 
attracting increasing attention in schizophrenia research. 
Two studies have experimentally validated five schizophrenia- 
associated genes as targets of miR-137 (Kim et al, 2012; 
Kwon et al, 2013). Through analysis of expression profiles 
in normal human tissues, we found that MIR137 is 
specifically expressed in the brain (Liang et al, 2007) 
(Supplementary Figure 1). Moreover, recent studies have 
revealed that MIR137 rsl625579 is associated with 
impaired cognitive performance of schizophrenia patients 
(Cummings et al, 2013; Green et al, 2013), and accumulating 
biological evidence, including the observed role of miR-137 
in regulating adult neurogenesis (Szulwach et al, 2010) and 
neuronal maturation (Smrt et al, 2010), supports the pivotal 
role of MIR137 in the pathogenesis of schizophrenia. 
Deciphering the neural mechanisms underlying the associa- 
tion between the MIR137 variant and schizophrenia may 
help us to establish the pathway from microRNA to risk for 
schizophrenia. Previous evidence has supported that miR- 
137 is widely expressed throughout the dentate gyrus of the 
hippocampus, suggesting its important functions in both 
developing and mature neurons (Smrt et al, 2010). There- 
fore, our finding of the influence of MIR137 on the DLPFC- 
HF connectivity supported our speculation that miR-137 
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affects the development and plasticity of the hippocampal- 
cortical circuitry. Considering that dysconnectivity of the 
DLPFC-HF circuit is viewed as one of the typical inter- 
mediate phenotypes in schizophrenia patients and in 
individuals at high genetic risk of the disease (Meyer- 
Lindenberg et al, 2005; Esslinger et al, 2009; Rasetti et al, 
2011; Paulus et al, 2013, 2014), our findings provide a 
possible neural mechanism through which MIR137 may 
impact brain connectivity and the risk for schizophrenia. 
The existence of different DLPFC-HF connectivity in 
individuals homozygous for the MIR 137 risk allele suggests 
that individuals who are at a high risk for schizophrenia 
show disrupted interactions between the DLPFC and the 
hippocampus. In fact, previous studies have consistently 
suggested that the DLPFC and the hippocampus belong to 
task-positive network and task-negative network, respec- 
tively (Fox et al, 2005; Hampson et al, 2010). The intrinsic 
functional brain organization may reflect the resting brain 
shifts between two different modes of processing. 

Two recent studies attempted to measure the effect of 
MIR137 on brain structure (Lett et al, 2013) and brain 
function (Whalley et al, 2012). Based on neuroimaging 
measures of brain structure, Lett et al (2013) reported that 
schizophrenia patients with the MIR 137 risk genotype 
exhibited reduced white matter integrity, as well as smaller 
hippocampi and larger lateral ventricles, whereas the brain 
structure of patients with the protective allele displayed no 
difference from that of healthy control subjects. Based on 
fMRI data during a sentence completion task in people at 
high genetic risk of schizophrenia or bipolar disorder, 
Whalley et al (2012) reported that the variant of MIR137 
influences brain activation in the posterior right medial 
frontal gyrus. Here, we examined resting-state functional 
connectivity, which can be used to explore the overall 
organization of the functional brain network (van den 
Heuvel and Hulshoff Pol, 2010) and may be altered in 
patients with psychiatric disorders (Hulvershorn et al, 2011; 
Rosazza and Minati, 2011). These findings presented above 
suggest the putative neural mechanisms underlying the 
relationship between MIR 137 and the risk for schizo- 
phrenia. Our results from healthy Han Chinese subjects may 
provide new clues for understanding the pathological 
mechanisms of schizophrenia. 

Another interesting finding of the present study was that 
individual DLPFC-HF connectivity affected WM performance 
in a different way depending on the MIR137 genotypes. For 
MIR137 TG individuals who may confer lower risk for 
schizophrenia, individual DLPFC-HF resting-state connec- 
tivity was inversely correlated with individual WM perfor- 
mance. Such correlation suggests that the individuals with 
higher DLPFC-HF connectivity have worse WM perfor- 
mance. In fact, the connectivity between task-positive 
and -negative brain regions has been reported to be 
associated with WM performance (Kelly et al, 2008; 
Hampson et al, 2010; Unschuld et al, 2013), which is one 
of the most typical cognitive impairments in schizophrenia. 
The disappearance of the inverse relationship of connectiv- 
ity and WM capacity in MIR137 high-risk individuals may 
suggest a disrupted modulation of DLPFC-HF connectivity 
with cognitive performance in TT individuals. 

As the MIR137 regulating the ZNF804A (Kim et al, 2012) 
and ZNF804A variants has been reported to impact 
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DLPFC-HF functional connectivity (Esslinger et al y 2009, 
2011; Rasetti et al y 2011; Paulus et al y 2013), we further 
examined the impact of the ZNF804A SNP (rs 1344706) on 
DLPFC-HF connectivity and also examined whether the 
current findings regarding MIR 137 were affected by the 
genetic variants of ZNF804A SNP found in the present 
sample. All of the current results for MIR 137 remain 
unaffected by controlling for the ZNF804A genotype as a 
new covariant in the group comparison analysis. Moreover, 
we did not find any significant differences in the resting- 
state DLPFC-HF connectivity or WM performance among 
the three ZNF804A genotypes (one-way analysis of variance, 
corrected P>0.05). This result may be explained by the 
findings of a previous study, which reported that ZNF804A 
influenced DLPFC-HF connectivity during the WM task but 
not during the resting state (Esslinger et al y 2011). 

Our findings may provide new clues to understand the 
neural mechanisms of MIR 137 on brain function and risk 
for schizophrenia. However, several issues should be 
discussed here. First, task-independent and task-dependent 
DLPFC-HF functional connectivity have been studied 
previously (Meyer-Lindenberg et al y 2005; Esslinger et al y 
2009; Rasetti et al y 2011; Paulus et al y 2013, 2014), whereas 
different seed regions were examined. Either the DLPFC 
region, which shows strong activation during specific tasks, 
or the HF region, derived from specific brain atlases, has 
been selected as the seed region in these works. Our 
findings of the present study were based on voxel-wise 
DLPFC functional connectivity maps throughout the whole 
brain using the DLPFC as the seed region. We further 
investigated the influence of different ROI on the DLPFC- 
HF functional connectivity analyses by taking left HF as 
the seed region (Supplementary Text). Voxel-wise group 
comparisons conducted throughout the prefrontal cortex 
also indicated that individuals homozygous for the MIR137 
risk allele (TT) display significantly different HF connectiv- 
ity with the right DLPFC (peak voxel MNI coordinate: 
x = 39, y = 33, z = 30; T = 3.71; cluster size = 43) relative to 
other individuals (Supplementary Figure 2). However, how 
to obtain the best seed region is still a question for the 
region-based functional connectivity analyses. Second, the 
difference of the sample size between two MIR137 geno- 
types is relatively large (TT = 249 and TG = 41), although 
simulation-based comparisons were performed in the 
present study by randomly resampling 41 of the 249 TT 
subjects and comparing the random subsample with the 
same number of TG subjects. The validation of our findings 
in independent samples is desired for the future study. 
Third, in addition to the disturbed DLPFC-HF connectivity, 
which is the focus of our current study, the functional 
connectivity among several other regions was also found 
to be different between the two groups (Supplementary 
Figure 3 and Supplementary Table 2). Fourth, our current 
analysis indicated that the DLPFC-HF connectivity was 
negative. Such negative functional connectivity of the 
DLPFC-HF was also reported in previous studies, especially 
for resting-state connectivity (Fox et al, 2005; Meyer- 
Lindenberg et al, 2005; Zhou et al y 2008; Rasetti et al y 
2011). This negative connectivity might be explained as 
functionally competitive interactions among the neuronal 
activity of different brain regions (Fox et al y 2005; Zhou 
et al y 2008; Buckner et al y 2013). However, the sign and 



magnitude of functional connectivity must be interpreted 
cautiously because it is a relative measure that depends 
heavily on the preprocessing and denoizing steps of fMRI 
data before performing functional connectivity analysis 
(Buckner et al y 2013). Recently, the relationship between 
global signal regression during the preprocessing stage of 
resting-state fMRI data and negative connectivity is still 
under debate. Several studies have suggested that the 
anticorrelation is most likely caused by regressing out the 
global signal (Murphy et al y 2009; Weissenbacher et al y 2009; 
Anderson et al y 2011; Saad et al y 2012). Other studies argued 
that the estimates of functional connectivity by not 
including global effects as confounding factors might yield 
bias, which obscures the underlying relationship between 
different brain regions, and negative correlations between 
default mode network (task-negative network) and task- 
positive network could be confirmed by other methods 
(Chang and Glover, 2009; Carbonell et al y 2011, 2014). To 
test the potential influence of global signal regression on 
our findings, we reanalyzed our data without global signal 
regression. The results without global signal regression 
indicated that the main findings kept consistent with the 
results with global signal regression. There were signifi- 
cantly different functional connectivity of DLPFC-HF 
between two genotypes (P = 0.041). Moreover, the DLPFC- 
HF connectivity without global signal regression also could 
predict the WM performance in MIR137 TG individuals 
(R= -0.437, P = 0.003), but not in TT individuals 
(#=-0.030, P = 0.322) (Supplementary Figure 4). Here 
we aimed to discover the group differences of the resting- 
state functional connectivity between different genotypes 
and its relationship with WM performance by considering it 
as a relative measure. In this opinion, the main findings 
with and without global signal regression kept consistent in 
our data. Nevertheless, although it was not the focus of our 
current study, we would like to note that the negative 
connectivity should be interpreted cautiously since we have 
observed the influence of global signal regression on the 
sign of resting-state functional connectivity. Finally, 
although neuroimaging evidence can increase our under- 
standing of the association between MIR137 and schizo- 
phrenia, the effects of MIR137 rs 1625579 on the expression 
level and the molecular function of its target downstream 
genes, as well as the consequent changes in the neural 
system function, remain unclear. Only one recent study 
(Guella et al y 2013) based on findings obtained through 
post-mortem examination has reported that miR-137 
expression levels in the DLPFC are decreased in homo- 
zygous MIR137 rs 1625579 TT subjects compared with TG/ 
GG subjects. However, additional biological evidence is 
needed to establish a causal relationship between MIR137 
rs 1625579 and the risk of schizophrenia. 
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